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In the presented paper known (up to the beginning of 2008) Lie- and non-Lie exact so- 
lutions of different (1 -I- l)-dimensional diffusion-convection equations of form f{x)ut = 
{g{x)A{u)ux)x + h{x)B{u)ux are collected. 

1 Introduction 

This is a review paper, where we present a brief summary of known exact solutions of variable 
coefficient (1 -|- l)-dimensional diffusion-convection equations of form 



where / = /(x), g = g{x), h = h{x), A = A{u) and B = B{u) are arbitrary smooth functions of 
their variables, f{x)g{x)A{u)^Q. 

Our aim is not to give a physical interpretation of the solution of diffusion equations (that 
is too huge and cannot be reached in the scope of a short paper), but to list the already known 
exact solutions of equations from the class under consideration. However, in some cases we give 
a short discussion of the nature of the listed solutions. 

The majority of the listed solutions have been obtained by means of different symmetry 
methods, such as reduction with respect to Lie and non-Lie symmetries, separation of variables, 
equivalence transformations, etc. 

Let us note that the constant coefficient diffusion equations (/ = 5 = 1, i? = 0) are well 
investigated and some of exact solutions given below were summarized before in [26,48]. 

Our paper is organized as follows. First of all we adduce solutions of the linear heat equa- 
tion obtained by means of various symmetry methods. In Section 3 the linearizable Burgers, 
Fujita-Storm and Fokas-Yortsos equations are considered. Lie reduction of constant coefficient 
nonlinear diffusion equation {KB = 0, / = 5 = 1) is performed in Section 4. Solutions of constant 
coefficient diffusion equations with exponential nonlinearity are adduced in Section 5. Solutions 
of constant coefficient diffusion equations with power nonlinearity are presented in Section 6. 
The important particular case of such equations, namely, the fast diffusion equation, is studied 
in more detail in Section 8. Diffusion equations with other nonlinearities are briefly discussed 
in Section 9. The next considered case (Section 10) covers the nonlinear constant coefficient 
diffusion-convection equations {f = g = h = 1). In Section 11 we adduce a brief analysis of 
known solutions of n-dimensional radially symmetric nonlinear diffusion equations. In Section 12 
exact solutions of some variable coefficient diffusion-convection equations are collected. At last, 
in Sections 13 and 14 we present a detailed analysis of interesting variable coefficient equations 
having distinguished invariance properties. 

In the Appendix A we adduce the complete results of group classification of equations (1) 
with respect to the extended group (5~ of equivalence transformations (22). 

Below, if it is not indicated separately, a, Sj, A, a, b, c, q are arbitrary constants, e = ±1. 
For convenience we use double numeration T.N of classification cases and local equivalence 
transformations, where T denotes the number of table and N does the number of case (or 
transformation, or solution) in table T. The notion "equation T.N" is used for the equation of 
form (1) where the parameter-functions /, g, h, A, B take values from the corresponding case. 





2 Linear heat equation 



Systematical investigation of invariant solutions of different dilfusion equations was started by 
the case of linear heat equation [32,33,39,42,44,45] 



Ut = Ua 



(2) 



which is invariant with respect to the six dimensional symmetry algebra generated by the vector 
fields 

Qi = dx, Q2 = dt, Q3 = udu, Qi = 2tdt + xdx, 
Q5 = 2tda: - xudu, Qe = At^dt + 4txdx - {x^ + 2t)udu 

(For the moment we are ignoring the trivial infinite-dimensional subalgebras coming from the 
linearity of the heat equation and corresponding to the linear superposition principle). 

The most general solution obtainable from a given solution u = f{t,x) by group transforma- 
tions is of the form 
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where v{t,x) is an arbitrary solution to the linear heat equation [42]. Considering the higher- 
order symmetry generators, one can prove that ii u = fit, x) is a solution of the linear heat 
equation than 

u = 2tfx{t, x) + xf{t, x) and 

u = t^ft{t, x) + txfxit, x) + J {2t + x^)f(t, x) 

are also solutions of the same equation [4] . 

All possible inequivalent (with respect to inner automorphisms) one-dimensional subalge- 
bras of the given algebra are exhausted by the ones listed in Table 1 [68] together with the 
corresponding ansatzes and the reduced ODEs. 

Table 1. Reduced ODEs for linear heat equation (2) 
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Thus we have the following solutions of (2): 
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« = cie" * cosh(Q!X -|- C2), u = cix + C2, u = cie " * cos(ax -|- C2), 
where U{b,z), V{b,z), W{c,z) are parabolic cylinder functions [1]. 



Multiplicative separation of variables leads to the solutions invariant with respect to {Q2 + 
aQz) ■ The additive separation of variables yields an exact solution of form 

u = cix^ + C2X + 2cit. 

The known Q-conditional symmetry operators and corresponding reductions are adduced in 
Table 2. [21] 

Table 2. Nonclassical reductions for linear heat equation (2) 



N 


Operator Q 


Ansatz 


Reduced equation 


1 


—xdt + dx 




V3" = 


2 


— Xdt + 9a: + X^du 




ifi" = -3 


3 


x^dt — Sxdx — 2>udu 


u = xtp (^t+ 


V3" = 


4 


x^dt — 2>xdx — (3w + x^)du 


U = Xip(t+^'^+^ 


^" = -15 
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xdx + udu 


u = xip(t) 


if' =0 


6 


coth xdx + udu 


u = (fi{t) cosh X 


if' -ip = ^ 
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— cot xdx + udu 


u = ip{t) cos X 


+ = 
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dx - udu - ii^du 


u = {2t- x)e-''ip{t) 


V?' - = 
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dx - ^-2{t + U)du 


u = ~t- ^[x + ip{t)]'^ 


93' =0 


10 


(t + 4) ft - xdx 


U = (fi(tX+ 


(^" = 



If function /(i, x) is an arbitrary solution of the linear heat equation and u is the general 
integral of the ODE j^dt + jdx = 0, then u satisfies the linear heat equation. This statement 
can be considered as another algorithm of generating solutions of the linear heat equation [21]. 
Indeed, even starting from a rather trivial solution of the heat equation u = 1 one gets the chain 
of quite interesting solutions 
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and among them the solutions 
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called often the heat polynomials [69]. 
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3 Linearizable equations 

Class (1) contains three equations, namely Burgers equation 

Ut = Uxx + 2wWa,, 

Fujita-Storm equation 
and Fokas-Yortsos equation 

Ut — {u~^Ux)x + U~^Ux, 



(4) 
(5) 
(6) 



3 



that are linearizable by the potential equivalence hodograph transformation and additional local 
equivalence transformations [10, 14, 18, 19, 25, 38, 50, 63, 64] : 



Uf = (u '^Ux)x + U "^Ua. 



t = i,x = v,u = u ^,v = x 



Ul = Uxx + 2uUx 



Vx = U, Vf = Ux+U 



r,2 



t = t, x = e^,u = e ^u, V = V 



u = —, Vt 

V 



t = t, X = X, u = ue", V = 













Ul = {U '^Ux)x 
Vx = U, Vf = U~'^Ux 


t = t, X = V,U = U ^, V = X 


Vx =U, Vt = Ux 





Therefore, applying the above transformations to the well-known solutions of the linear heat 

equation one can easily construct solutions of the linearizable equations. 

Thus, e.g., the fundamental (source) solution u = (47ri)~^/^e"^' Z^"^*) and dipole solution 
u = — ((47rt)~^/^e~^ ^^^^^)x are mapped into the separable and self-similar solutions of the Fujita- 
Storm equation [48] 



u = {Anty/^e"'' where x = tt'^/^ f e'v' 

Jo 



dy, 



-1/2 



correspondingly. 

Other solutions of the linear heat equation presented in the previous section yield the following 
explicite exact solutions of the Pujita-Storm equation [30]: 



u = c, u 
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1 
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In [66] the following formula for deriving exact solutions of the Fujita-Storm equation is 
derived: if u{t, x) is a solution of the Fujita-Storm equation (5) then 



V = u + 



2tUt + U + XUx 
-t{Ux — |tt^) — f W 



is also solution of the same equation. 

Similarly [67], if u(t,x) is a solution of the and the Fokas-Yortsos equation (6) then 



u 



V = 



u — u. 



X 
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is also solution of the same equation. In [67] a new exact solution of the Fokas-Yortsos equation 
is presented: 

1 t^(e^*+^) 



2 1 + VF(e4*+^) ' 

where W{x) is the Lambert W function determined as W{x)e^^^'> = x. Let us note that other 
exact solutions of the Fokas-Yortsos equation (6) can be easily recovered from the solution set 
of the Fujita-Storm equation (5) by means of application of the local transformation of variables 
shown in the above scheme. 

The same tricks can be used for obtaining exact solutions of the remaining linearizable equa- 
tions. However, since the adduced transformations are nonlocal, sometimes it could be easier 
to search directly for solutions of the nonlinear equations. Thus, e.g., one can easily find Lie 
solutions of the Burgers equation: 

ci — X 1 2x + ci 

u = - , u = £-\ , u 



2t + C2 x + 2et + c x2 + cix + 2t + C2 ' 

6 (x^ + 2t + c{) ci a e"^ (^^-^et) _ 

UL = U '=■ XL — S ~\~ 

2x3 + 12ta: + 6cit + C2 ' i + c2e-^?*-^i^ ' 2 e'^i(=^-2^*) + C2 ' 

Acos(Ax + ci 



u = —£ + ci tanh(ci(x — 2et) + C2), u = —-5 ( 2tanh ^ — Ax — C2 

2(A^t + ci) \ XH + ci 



-£ — ci tan(ci(x — 2rf) + C2), 



C2e-^^* + sin(Ax + ci) ' 
1 



ci / x + cs W X + C3 

u = — = exp — — cierr — , + C4 

V7r(t + C2) V 4(t + C2);V 2Vt+^ 

cos 2xe~3* + ci sin x cie~*(cos x + sin x) + 026*+^ 



— cosxsinxe 3* + ci cos x + C2e* ' c\e *(cos x — sinx) — C2e*+^ + 6 

Here erfz = e~^^dS, is the error function also called the probability integral. The last two 

solutions were found in [54]. 

Solutions of the Fujita-Storm and Fokas-Yortsos equations can be singled out form the 
solutions adduced in Sections 6 and 10 taking /i = —2, u = —2. (Note that all these solutions 
can be also reconstructed from ones of the linear heat equation by means of potential equivalence 
transformations . ) 

If u{t, x) is a solution of the Burgers equation (4) then 

v = u-\ u = u + 



Uj; + ' Ux + + u 

are also solutions of the same equation [66]. 

Potential equivalence transformations were used to obtain solutions of some boundary-value 
problems adduced in [10,41,63,64] 

4 Nonlinear diffusion equations. General case 

Consider now the class of nonlinear diffusion equations 

Ut = {A{u)Ua,)a., (7) 

where Au 7^ 0. Lie symmetries of this class have been studied in [44]. The Lie symmetry algebra 
of equation from class (7) with arbitrary value of parameter-function A{u) is three-dimensional 



5 



and spanned by 



Qi = dt, Q2 = dx, Q3 = "^tdt + xd^. 

Taking into account discrete symmetry transformations of changing sings of independent 
variables one can formulate the following statement. If ti = f{t, x) is a solution of equation (7), 
then u = /(ef i + £25 + £3) is also solution of the same equation. 

All possible inequivalent (with respect to inner automorphisms) one-dimensional subalgebras 
of the given symmetry algebra,the corresponding ansatzes and the reduced ODEs are exhausted 
by the ones listed in Table 3. 

Table 3. Reduced ODEs for nonlinear diffusion equations (7), ^ 0. 
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The first three equations can be easily integrated for all values of A{u). Solutions of the last 
equation are known for many functions A{u) (see Section 9). 

Consider now in more detail equations with wider symmetry algebras. Up to the group of 
equivalence transformations 

t = eit + 64, X = e2X + £5, u = £su + £e, A = £'^^£\A 

there exist three inequivalent cases of extensions of Lie symmetry algebra [44]: a = e", a = u^, 
H + -4/3 and a = u'^l'^ . 

5 Nonlinear diffusion equations. Exponential nonlinearity 

We start from the equation with exponential nonlinearity 

ut = (e"wa,)a„ (8) 
having the four-dimensional Lie algebra spanned by the operators 

Qi = dt, Q2 = tdt - du, Q3 = dx, Qi = xdx + 25„. 

The only non-zero commutators of these operators are [Qi, Q2\ = Qi and [Q3, Qa\ = Qs- There- 
fore A"^'^^ is a realization of the algebra 2^2. 1 [40]. All the possible inequivalent (with respect 
to inner automorphisms) one-dimensional subalgebras of 2^2. 1 [46] are exhausted by the ones 
listed in Table 4. 



Table 4. Reduced ODEs for (8). a 7^ 0, e = ±1, S = signt. 
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Ansatz u = 




Reduced ODE 
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(Qs) 
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(Q4) 


ifiijj) + 21n|a;| 
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5 


(Qi + eQa) 
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X — et 
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(Q2 +£(33) 


ifiiu}) — In \t\ 
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(e^)" = -6{e<p' + 1) 
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Solving the equations 4.1-4.5 we have the following solutions of (8): 

, - ,-2 

u = In Icix + Co|, II = In 



2t 



+ 



-dip = X — et + cq. 



u = ipix — et) where 

■ ci-eip 

Equation (8) admits an additive separation of variable that leads to the solution invariant 
with respect to scale transformation. 

liu = f{t, x) is a solution of equation (8), then 

U = f{eit + £3, £2^; + £4) - £1 + 2^2 

is also solution of the same equation. 



6 Nonlinear diffusion equations. Power nonlinearities 

Another case of equations admitting extension of the Lie symmetry algebra is the one having 
power nonlinearity 

wt = (|nrw,),. (9) 

As in the previous cases the invariance algebra of (9) 

A"^^ = (Qi = du Q2 = tdt - n'^udu, Qs = d^, Qi = xd^ + 2n-^udu) 

is a realization of the algebra 2^2.1 • 

If = f{t,x) is a solution of equation (9), then 

u = £^^£2/(^1 ^ + ^3' ^2^ + ^4) 

is also solution of the same equation. 

The result of reduction (9) under inequivalent subalgebras of is written down in Table 5. 

Table 5. Reduced ODEs for (9). ^ / a 7^ 0, £ = ±1, 5 = signt. 



N 


Subalgebra 


Ansatz u — 
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(Qi) 
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X — et 
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For some of the reduced equations the general solutions are known. For other ones we 
succeeded to find only particular solutions. These solutions are following: 

1 



u = \cix + co\f+'^ , u = (cq — £ijl{x — et))i^ , u 



H [x + co) 



H + 2 



2t 



+ Cl\t\ M+2 



u 



11 + 2 2t 
u = (f{x — et) where 



+ ci(a; + co)''+i|i| ^(m+i)^ 



J ci- 



eip 



-dip = X — et -\- cq. 



(10) 
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Equation (9) admits multiplicative separation of variables. Namely, for all values of /x one 
can find solution in form of the product of two functions of different arguments: 

u{t,x) = {ht+bor^^''f{x), (11) 

where the function / = f{x) is given implicitly 



/ 



= ±x + Co, A 



These could be found also in [2,5,6,15,24,26,48,52,58,72,74]. The most studied cases 
are equations with fi = ±1, —2, —4/3, —3/2. Below we adduce their exact solutions that are 
inequivalent to (10) with respect to the Lie symmetry transformations. 

Equation with the singular value of the parameter n = —1 called often the fast diffusion 
equation, is distinguished by the reduction procedure. Lie invariant solutions of it will be adduced 
in separate section together with non-Lie solutions obtained from invariance of the fast diffusion 
equation with respect to the nonclassical potential reduction operators. 

Fujita-Storm equation (5) is linearizable and has been considered in a separate section. 

Equation 



Ut 



{u-^/^u^)^ (12) 



admits the five-dimensional Lie symmetry algebra generated by 

3 

Qi = dt, Q2 = tdt - iT^udu, Qs = dx, Q4 = xd^ + -udu, Q5 = -x^d^ + 3xu5„. 
An optimal system of one-dimensional subalgebras of this algebra is 

(Q1 + Q4), {aQz + QA), (Q5), (<32 + (35), {Qs + Qb)- 

One can easily construct the corresponding ansatzes and the reduced ODEs. However, to the best 
of our knowledge all the found solutions of these equations are equivalent to (10) with particular 
value of parameter /x = —4/3. Solutions of (12) invariant with respect to dilatation operators 
can be found also in [35]. Besides the already adduced Lie invariant solutions, equation (12) has 
functional separated solution [22, 57] 

u = {<p4it)x^ + fz{t)x^ + ^2{t)x^ + Mt)x + <^o(^))"^/^ 
where the functions cpi = (pi {t) are determined by the system of ordinary differential equations 

3 3 
</^o = -^V'l + 2<^o<^2, y^'i = -^1^2 + Qfo^a, y^'2 = -^2 + 2'^W3 + '^'^^o'Pa, 

3 

f's = -'P2'P3 + Q'Pl'PA, 9?4 = --9?! + 2(^2¥'4- 

The general form of exact solutions of (12) obtained from the known ones u = f{t,x) with 
action of group transformations is 

-12 / -4/3 \ 

£1 „ / -4/3 , . £2 ^ 



{£5X -I- 1)3 \^ ^ e5X + l ^ 

Equation 

Ut — (u'^/^Ux^^, 
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admits also the functional separation of variables. The corresponding exact solution is 

u = (Scix^ + f2{t)x^ + Mt)x + Mt)r^^\ 

Here 

f2it) = 3j^it)dt + 3C2, flit) = -(/ ^it)dt + C2)2 + ^^(t), 

Ci ZCi 

fo{t) = ^(/ f{t)dt + C2)3 + V{t)dt + C2) + ^^V'{t), 

where the function (p = (p(t) is defined implicitly by / (03 — 8(p^)~^^'^d(p = ±i + C4 . 
T.K. Amerov [3] and J.R. King [36] suggested to look for solutions of the equation 

in the form u = {ip^{x)t + where the functions ip^{x) and (^^(a;) satisfy the system of 

ODEs tp].^ = (ip^)'^, iPxx = 'P'^f^- ^ particular solution of this system is 

16 Ci C2 

<^ = — ' = — + 

X 3C X 

7 Porous medium equation ut = {uUx)x 

Another important subclass of diffusion equations is a special case of equation (9) with /x = 1 

Ut — {uux)x, (13) 

called also porous medium equation. It first exact solution has been obtained by Boussinesq [13]. 
He was looking for a solution in a separated form u(t, x) = X{x)T{t) satisfying conditions 

u{t,0) = 0, U^\x=L = 0- 
Thus constructed solution reads as 

Horn 

l + {3b^Ho/2L^)t' 

where Hq is a constant, ^ = x/L and the function F = F(^) is defined implicitly 



1 Xd\ XdX _l„/2 In 

bJo TT^' Jo 7r^"3 U'2J- 



^ = 1 

The next exact solution of the porous medium equation (13) was found much later by Baren- 
blatt [6] and written in the present form by Sokolov [59] (the instant source solution): 

„ = l((9t)2/3_^2-)^ < X < (9i)V3 = 

It is easy to become convinced of the fact that both Boussinesq and Barenblatt solutions cor- 
respond to the Lie symmetry of the porous medium equation (13) [73]. Indeed, the Boussinesq 
solution is a particular case of the ansatz 

u={l + at) (p{x), a= , 



9 



that is invariant with respect to the one-parameter Lie symmetry group generated by Q = 
{l+at)dt — udu- Barcnblatt solution is invariant with respect to the one-parameter Lie symmetry 
group generated hy Q = Stdt + xd^ — ud^- 

Besides (10), its exact solutions in parametric form are known [48]: 

X={6t + Ci)^ + C2f + C3, U= -{6t + Ci)C^ - 2C2^^ 

X = tf{u;) + g{uj), u = tf{uj)+g'{uj), 
where the functions / = /(w) and g = g{uj)ave determined by the system of ODEs: 

if? - ff" = /'", f'g' - fg" = 9"'. 

It is obvious, that the second equation has two linearly independent particular solutions g = 1 
and g = f ■ The general solution of these equations can be represented in form 

5 = ci + C2/ + C3(//Vrfc^-//«, / = /H, ^=(^e-//<^. 

It is not difficult to verify, that it has the following particular solutions 

/=^— and / = cie^^'^. 

U) + Cl 

One can see, that the first solution leads to the previously given implicit solution. 

8 Fast diffusion equation ut = iu~^Ux)x 

All invariant solutions of fast diffusion equation 

Ut = (u-^Ur,)^ , (14) 

which were earlier constructed in closed forms with the classical Lie method, were collected 
e.g. in [48-51]. A complete list of Gi-inequivalent solutions of such type is exhausted by the 
following ones: 

(15) 

4)^ = ^-—^' 5)u= y—, 6)u = — -g— , 7) u 



x"^ + et"^ ' cos^ X ' cosh^ x ' sinh^ x 

The below arrows denote the possible transformations of solutions (15) to each other by means 
of the potential hodograph transformation (16) up to translations with respect to x [50]: 

O l)e=0 ; 1)2=1 < > l)e=-l, a;+t<0 ; O 1)£=-1, x+t>0 ; 2) < > 3)^^=0, a;>t ; 

O 4)^=0 ; 5) i — > 4)^=4 ; 6) < — > ^)e=-A, \x\<2\t\ ; 7) < > ^)e=-A, \x\>2\t\ ■ 

The sixth connection can be found also in [20,53]. If /x 7^ solution 3) from list (15) is mapped 
by (16) to the solution 

u = Mioj) - t + nte~^^'^\ uj = x-\n\tl 

which is invariant with respect to the algebra {tdt + dx + udy)- Here i) is the function determined 
implicitly by the formula J"(i? — 1 + iie~'^)~^d'd = u. 
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Some non-Lie exact solutions of (14) were obtained in [23,55,56]. Thus, P. Rosenau [56] 
found that potential equation vt = Vx~^Vxx corresponding to (14) admits, in addition to the 
usual variable separation v = T(t)X{x), the additive one v = Y{x + At) + Z{x — Xt) which 
is a potential additive variable separation for the fast diffusion equation (14). (The classical 
multiplicative separation of variables is given by (11) with /x = — 1.) To construct nonclassical 
solutions of (14), C. Qu [55] made use of generalized conditional symmetry method, looking for 
the conditional symmetry operators in the special form Q = {uxx + H {u)ux^ + F {u)ux + G{u))du- 
M.L. Gandarias [23] investigated some families of usual and potential nonclassical symmetries 
of (9). In particular, using an ansatz for the coefficient 77, she found non-trivial reduction 
operators in the so-called "no-go" case when the coefficient of dt vanishes, i.e. operators can be 
reduced to the form Q = dx + r]{t,x,u)du- These solutions and the ones similar to them were 
represented uniformly over the complex field as compositions of two simple waves which move 
with the same "velocities" in opposite directions in [51]. Using such representation the following 
solutions of fast diffusion equation (14) were obtained [51]: 

, , , , 2sin2t 

1 ) u = cot[x — t) — coiyx + 1) 



2') u = coth(a; - t) - coth(x + t) = 
3') u = coth(a; — t) — tanh(a; -|- t) - 
4') u = tanh(a; — t) — tanh(x + t) ■ 
5') u = cot(ix + t) — cot (ire — t) = 
6') u = i cot{x + it) — i cot{x — it) 



cos 2t — cos 2x ' 
2 sinh 2t 



cosh 2x — cosh 2t ' 

2cosh2i 
sinh 2x — sinh 2t ' 

2 t,iuh 2L 
cosh 2x + cosh 2t ' 
2 sin 2t 
cosh 2x — cos 2t ' 
2 sinh 2t 



cosh 2t — cos 2.x 

Transformation (16) acts on the set of solutions l')~6') in the following way [51]: 

l')cos2t<cos2a; ■* ^ ^')\t-^t+Tr/2,x-*x/2,v^2v'i l')cos 2t>cos 2x ^ ^')\x->-x/2,v^2v-it ! 
2')|x|<|t| ^')\x-*x/2,v-*2v'i O 2')\j.\^\t\\j._^^/2,v-*2v 'i 

O ^')x<t\x^x/2,v^2vj ^')x>t ■* ^')x>t\x^-x/2,v^-2v'i O G')\x^x/2,v^2v 

These actions can be interpreted in terms of actions of transformation (16) on the nonclassical 
symmetry operators which correspond to solutions l')-6'). 

In [56] P. Rosenau considered additive separation of variables for the potential fast diffusion 
equation (14) and constructed solution 4'). Using the generalized conditional symmetry method, 
C. Qu [55] found solutions which can be written in forms 1') and 6'). After rectifying compu- 
tations in two cases from [55], one can find also solutions 2') and 5'). Solutions 1'), 3') and 4') 
were obtained in [23]. The remaining solutions from the above list were found in [51]. 

One of techniques which can be applied for finding the above solutions is reduction by con- 
ditional symmetry operators of the form Q = dx + [v^it, x)u + rf{t, x))udu (see [23] for details). 
All reductions performed with reduction operators of such type result in solutions which are 
equivalent to the listed Lie solutions l)-7) or solutions l')-6'). 
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9 Nonlinear diffusion equations. Other nonlinearities 



Known exact solutions of the redueed equation 3.4 corresponding to (Q3 = 2tdt + xdx) (which 
are self-similar solutions of (7)) are adduced in Table 6 [48]. 

Table 6. Self-similar solutions for nonlinear diffusion equations (7), oj = xj^/t 



N 



10 



11 



A{u) 



n „ n 
—11 1 

2 2{n + l) 



2{n + l) 



((1-M) 



n-l 



(1-nf") 



2(1 -n) 



1 . 2 TTW 

1: Sin — — 

2 2 



— sin TTit {-Ku + sin -ku) 
8 



16 



sin ITU (5 + cos iru) 



_cos- (^COS- + — -1 



ti arccos u + 1 1 
2v/l - u2 2 



TT — 2(1 — u) arcsin(l — u) T" 
4V2m - m3 2_ 



u arcsm u i o 



Solution u) = u!{u) 



{l-uT 



- 1 



1-sm- 



arcsin(l — u) 



VI 



- \nu 



Conditions 



n > 



n > 



< n < 1 



Let us give some more examples of travelling wave solutions [26,48,58]: 
Ut = (sinh^ uUx)x, u = arccosh^ ^ ~^ \ , 
Ut — (cosh^ uUx)x, u = arcsinh^ '^2t ' 

26ci 



Ut {{u^f^ + buf^)ux)x, u = 
Ut = ((e^" -I- bue^)ux)x, u = In 



ci 

(Cl + l)(c2-i) 
±X + Cl 



{X + Cl) - 



a(ci 1) 



Vc2 - 2t 

Ut — {ue'^Ux)x, u = hi{cix + eft + cq). 
For equation with logarithmical nonlinearity 

Ut = (lnuux)x 

travelling wave and self-similar solutions are known: 



u 



exp I lb a/ 2cix + 2cft + C2] , u = exp 



±X -|- Cl 

C2 -2t 



1 



1/m 



A number of exact solution for equations of class (7) were constructed with nonlocal (quasilo- 
cal or potential) symmetries [2,11,12,60-62]. 

Thus, e.g., reductions with respect to the optimal system of subalgebras of Lie algebra of 
potential/quasilocal symmetries of equation 

Ut — ((1 + U^)~^Ux))x 



12 



give rise to exact solutions of form 

u = ce*--(l - c2e2(*-))-V2^ „ = -2t- x'r^/\ 

u = tan((^(a;) + arctan(A(a;)) + £t), where uj = x'^ + v^, u = tan((^(a;) + st), 

9 9 

u = xtan((^(a;) + arctan(A(u;)) + et), where oj = , v = xtan + "2^*J • 

Here (p{u)) and A(ci;) arc arbitrary solutions of the system (p' = lu~^X/2, X' = (1 + A^)(e — aj~^A)/2. 
The list of known exact solutions of this equation (called often the Pujita's equation) involves 
also the following ones [2,16,26,48]: 

u = tan(cix + C2), u = ±a;(ci — 2t — x^)"^/^, 

u = — =, e(x + et) + Co = In , + ci arctan u 

^1 _ e2(t-^) ^ ^ 4 + 1 \ ^/WTl 

u = —^^=^= , where v = . = or 

vl-v^ V4ciC2 + C3e-2A2t 

ci sin Ax + C2 cos Ax 



q + C2 + cae^-^ * 

siiili x ± cosli sill x 

u —— u —— u —— 

V- cosh^ X - e-2* ' V- sinh^ x + e-^* ' Vcos^xi e^* 

We adduce also some exact solutions of another equation with Fujita's type nonlinearity 

Ut — ((1 - U^)~^Ux))x, 

namely [16]: 

X 3: 

u = c, u = tanhx, u = = (i > 0), u 



sinhx sinhx , , it cosh x 
u = , u = {t > 0), u= {t < 0), 

V cosh^ X + e~2* v cosh^ x — e"^* v sinh^ x + e^^* 



smx 



V— cos^ X + e^* 

The third similarity solution converges (pointwise) to a step function as t ^ O"*", and to zero 
as t — > 00. The fourth solutions are bounded travelling waves. The fifth and sixth solutions 
converge to the time independent second solution as t ^ 00. The sixth solution converge to a 
step function as i — 0+. The seventh solution converges to the values ±1 as f — 0— . The last 
solution converges to a square wave as t ^ 0+, and to zero as t ^ 00. 

One more example of solution obtained with application of potential symmetry is u = 
tan((/?(u;) + arctan(2t(;(/9') — A~-^ \nt) for the equation 

Ut = ((1 + n2)-ie^"'^^*"""w^))^. 

Here 

no l,c — w lb' lb 2 

u = x+v^ ^ = -ln^-arctan^(a;), + 2^ + = 

All the potential symmetries of equations from class (17) can be obtained from Lie symme- 
tries of (17) by means of prolongation to the potential variable v and application of potential 
equivalence transformations [50] 

i=t, x = x + £v, u= — ^ — , v = v, A = {l + eu)'^A 

1 + eu 
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and hodograph transformation 



i = t, X = V, u = u ^, V = X, A = v?A, (16) 

where = u, vt = Aitj - 

Therefore, these transformations can be used for obtaining potentially invariant exact so- 
lutions from the Lie ones. The complete list of nonlinear constant coefficient diffusion and 
diffusion-convection equations having potential symmetries together with the transformations 
mapping them to the equations with power and exponential nonlinearities can be found in [50] . 

In [65] non-point nonclassical symmetry operators are used to obtain exact solutions of evo- 
lution equations. In particular, it is shown that equation 

admits an exact solution in implicit form u = [z^ + c)/ (2t) , where 
z = — , if c> 0, and 

^ = , if c<0. 



1 z^+c I , ^-V^ I ci 



A number of authors considered additive separated solutions of diffusion equations (7), i.e., 
solutions of form 

u{t,x) = ip{t) + tp{x). 

Usually such solutions are Lie invariant. They were adduced in previous sections. So, here we 
adduce only list of equations admitting such kind of separation of variables. Namely, a diffusion 
equation (7) admits separation of variables if and only if it is C-equivalent to equation with 
the diffusion coefficient being of the following functions [16]: 

A = \u\'', ^ = e", A = {u'^±l)-\ 

A = z(u)e^(") , u = J^' s-^/^e-'/^ds, z > 0, 

A = e^^(") cosh z{u), a / ±1, u = [ cosh^^/^ se'^'^/'^ds, -oo < z < oo 

Jo 

^ = e^^(")sinhz(u), aj^±l, u = ^ sinh-^/^ se'^'/^ds, z > 0, 

A = 6"'^''^ cos z{u), u = [ cos-^^^ se-^'/^ds, -7r/2 < z < 7r/2. 

Jo 

10 Constant coefficient diffusion— convection equations 

Lie symmetries of the constant coefficient diffusion-convection equation 

Ut = {A{u)ux)x + B{u)ux, (17) 

B ^ and corresponding Lie reductions were considered by a number of authors, see for ex- 
ample, [17,32,33,43,71]. However the complete group classification of class (17) was presented 
only recently in [49] . 
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Any equivalence transformation of class (17) has the form: 

t = te|e5 + ei, X = xe4 + ert + 82, u = ue6 + e3, 
A = Ae^^, B = Bel^el^ - £7, 

where ei, . . . ,£7 are arbitrary constants, e^e^e^ 7^ 0. Note, that all the equations with B{u) = 
const are reducible to diffusion equation (7). Besides such 'trivial' cases, Burgers equation and 
Fokas-Yortsos equation, only few Lie invariant solutions for equations with non-zero convectivity 
are found. Thus, e.g., in [71] scale- invariant solution 

xu-^"'"^ - (1 + 2//i) ln(x«-''/2 - 1 - 2/ii) = cit + Co 

in implicit form is found for the equation with power nonlinearities 

was found. 

M. Edwards [17] investigated Lie symmetries of (17) and constructed optimal subalgebras of 
the symmetry algebras for some of equations from the class. Here we supplement her results 
and adduce the complete list of Lie reductions of equations from class (17). (The linearizable 
Fokas-Yortsos and Burgers equations have been considered separately in Section 3.) 



Table 7. Reduced ODEs for ut = (e^^^u^)^ + e"«^ (/x ^ 0) 



Subalgebra 


Ansatz u = 




Reduced ODE 


(dt) 

{dt + ed^) 

{{fxt -2t + s)dt 

+{fj, - l)xd^ + du) in =/- 2) 

{edt + xd^ + 9u) (a* = 2) 


9?(a;) 
V3(a;) 

ln((M - 2)i + e) 
fx~2 

•fiioj) + - 


t 
X 

X — st 
x{{n-2)t + e)^ 

xe-'l' 


V?' = 

(e^'^V'')' + e'^'P' = 

- + 1 = (e'"^(p')' + e-'^f' 



If u = /(t, x) is an exact solution of equation ut = (e^^Wa;)^ + e^u^-, then 
u = /(e('^-2)^ii + £3, e^'^-i^^i+^^x + £4) + £1 + £2 



is also solution of the same equation. 



Table 8. Reduced ODEs for ut = {e^Ux)x + uu. 



Subalgebra 


Ansatz u = 


UJ 


Reduced ODE 


{dt) 

(dt+ed,) 

{{t + e)dt + [x - t)d^ + du) 


<p{u)) 
ip{Lj) + ln\t + e\ 


t 

x 

X — et 


ifi' = 

ie'^ip'Y + = 

1 + <p'(i -uj) = {e^ip'Y + w' 



If u = f{t,x) is an exact solution of equation ut = {e^Ux)x + uux, then 
" = /(e^^* + £2,e^^{x — £\t + £3)) + £1 



is also solution of the same equation. 
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Table 9. Reduced ODEs for ut = {u^Ux)x + u'^Ux 



Subalgebra 


Ansatz u = 




Reduced ODE 


{da:) 

(dt) 

{dt + ed^) 

{{lit - 2vt + e)dt 

+ {p — u)xdx + udu), (m / 2;/) 

(edt + vxdx + udu), (p = 2j/) 




t 

X 

X — et 
x{{fj,- 2v)t + e)'^=^ 

xe-^'l' 


(/J = 

(lp^^)' + (^'V = 



If u = f{t, x) is an exact solution of equation ut = {u^Ux)x + u'^Ux, then 

u = e'^+''f{e^^-^''^'H + e3,e^^-''^''+''''x + £4) 

is also solution of the same equation. 

For equations with u = jj, the generalized travelling wave solution is known: 

_ (c2-x ln|t + ci|y/'' 

If = 2, = 1 then the generalized travelling wave solution in implicit form is the following: 

u^du 



-V? — 2eu + ci 



X — et + C2- 



For such values /v, and u the degenerate solution linear in x has the form u = T(i)(x + ci), where 
function r is given in implicit form 



+ 2 In 

r 



2r + l 



= t + C2. 



Table 10. Reduced ODEs for ut = {u^^Ux)x + Intttt^ 



Subalgebra 


Ansatz u — 




Reduced ODE 


{di) 

{d. + edt) 

{(fit + e)dt + {ixx — t)dx + udu) 


ifiioj) 


t 

X 

X — et 

u^x + e 1,1 . 


if' =Q 

(ipi^ip)' -\-\a.ipip' = Q 
-eip' = {tp'^tp)' + Inifiifi' 

—tp — jMjjip' = {(p'^ifiY + \n(p(p' 
A* 



If u = f{t,x) is an exact solution of equation ut — {u'^Ux)x + InwUa,, then 

u = e^^f{e''^H + £2, e''^' {x - eit + £3)) 

is also solution of the same equation. 

In particular, if = 0, we obtain ut = Wa;a; + ln uUx- This equation has two known travelling 
wave solutions (usual and generalized ones): 



u{t, x) = exp(cie ^^'''^^ + 1 - C2), u{t, x) = exp 



ci — X ^ In |f + C2I 



t + C2 t + C2 

Generalized travelling wave solution are known for the following equations: 



\ 1/m 



Ut = {{u^^ + bu'^)ux)x + ui^Ux, u = (^x(p{t) + ci(p{t) + ^(p{t) J (p{t)d?j 
Ut — ((ae^" + be'^)ux)x + e'^Ux, u = ln{x(p{t) + ci(p{t) + b(p{t) J (p{t)dt) , 
where ip = ip{t) is determined by the equation ip' = aip^ + ip^. 
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11 n-dimensional radially symmetric nonlinear 
diffusion equations 

Class (1) contains a subclass of physically important n-dimensional radially symmetric nonlinear 
diffusion equations. Preserving tfie common terminology we use for them the notation 



Equations of the form (19), especially for power nonlinearity 



(19) 



(20) 



have a large number of applications, for both /x > (slow diffusion) and /x < (fast diffusion). 

To the best of our knowledge, first similar solutions of (19) were instantaneous source-type 
solutions [6,47] namely, solutions of form 

if /X 7^ -2/N and 

where A is an arbitrary constant, in case fj, = —2/N. Substituting these values into (20) yields 
(assuming t > in the first three cases): 



M > : 
/X = : 

> |tx > -2/n 
fj, = —2/n : 

H < -2/n : 



u= < 



u = 



0, r>atV^^'*+2 



U = t-n/ [nn-\-2) 



2(/in+2) 



(^2 _^^2^-2/(;un+2)) 



u = e 



u = 



-Xnt 



-n/2 



(_i)-n/(/*n+2) 

0, 



^(a^ + r^(-t)-V(Mn+2)) 



1/m 



t < 
t > 



The following instantaneous source-type solutions for equations (20) were obtained by 
King [34]: 

Case /X = -1, n 7^ -2: / = (/xn + 2)exp ( - ^^j/J^exp ( - 4^)^^. 

Case Ai = -1, n = 1: / = l/(/?e-"^ - a-'^{iin + 2)~i(l + a^)). 

Case /X = -1, n = 2: / = l/(/3r + C^(^n + - ")"^)> if « 7^ 2 and 

/ = l/{pe + eil^n + 2)-i InC, if a = 2. 

Case /X = -1, n = 3: / = 2/(2/3e-°^/« + (^(a + + a~^e-°/«£;i(-a/C))(/xn + 2)-^), where 
£'1(2;) = e~^t~^dt is the exponential integral. 



Case /X = —2/n: / = ^ "5, where J (ra^i — ay^/" — \g^^'^/"-)dg = In^. 
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Case /i = —1/2: 



^l/2^(4-n)/2 

(//n + 2)1/2(4 _n) 



(4-n)/2 



l/2c(4-n)/2 



72^(4-^n)/2 



.(/xn + 2)1/2(4 -n) 



+ (l-/3)>^( 



72^(4-n)/2 



+ 2)1/2(4 -n) 

-2 



(/xn + 2)1/2(4- 



n 



where i/ = 2/(4 — n), J,y(^;) and i^j/(2;) are Bessel functions of the first and the second type 
corr esp ondingly. 

We now consider solutions [34] that generalize the one-dimensional dipole solutions given for 
/X > in [7]. The similarity variables are chosen to fix rudr in time (if the integral exists). 
In one dimension this corresponds to the centre of mass. The appropriate similarity solution 
to (20) then takes the form 

„^i-l/(M+l)/(^)^ ^ = ^i-l/(2(M+l))^ 

u = e-2^*/(6, ^ = re-^\ n = -1. 
For different values n and /x the following solutions are known [7, 34] : 



Case /X / 0, -2/n: / = ^(2-«)/(m+i) (^^^{(3 - 
Case 11 = 0: f = /^^^-ng-eV^ 

Case /x = -2/n: / = ( In i) 
Case n = -1/2: (here i/ = 2 - n/2) 

a \ /ai/2£N /ai/2£ ^ 



l/n 



ai/2e 



Case /X = -1/2, n = 3 (cr [70]): / = ^ tan^ (^^{^ - c)) or / = ^^tanh^ (^(^ - c)) . 



Case /X = —1/2, n = 1: / 



(2-0:1/2^ cot[al/2(^_c)/2])2 

Case /X = n/2 - 2: / = ^(2-»)/(;^+i)c,, where 

^n/2-2 2 



or / 



(2-aV2^coth[aV2(^_c)/2])2 ' 



5 



re 



n-2 



5 + (n-2)a (n-2)2 

12 Variable coefficient diffusion— convection equations 

To obtain invariant solutions of the variable coefficient diffusion-convection equations of form (1) 
two approaches were used. The first one is the direct finding of solutions invariant with respect to 
a subalgebra of the Lie invariance algebra, and the second one is reconstructing of new solutions 
from the known ones using equivalence transformations. 

As an example of implementation of the first approach we adduce some of the invariant 
solutions of equation [24,48]. 



\x\Put = {\u\f^u^)^. 



(21) 
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Namely, 



{p + 2){2 + n + p + np)J 

u 



\2 + p ) 



PI1 + P + lx + 2 



For equation (21) withp = — (/z + 2)/(/x + 1) solution obtained with multiplicative separation 
of variables is known: 

u = e-^i* (ci(// + i)2^M/(M+i)eCiMt + C2y^^ 

See also the next two sections for invariant solutions of essentially variable coefficient equa- 
tions. 

Another possible way of finding exact solutions is based on application of cqTiivalcncc trans- 
formations. The complete extended equivalence group of class (1) is formed by the trans- 
formations [27,28] 

t = 5it + 62, X = X{x), u = dsu + 6i, 

f = ^^f, g = eie^^X^ipg, h = eie^\h, 1 = ^2^ B = e^{B + e^A), (22) 
where 5j (j = 1,4) and {i = 1,4) are arbitrary constants, 5i5^eie2£^ 7^ 0, X is an arbitrary 

^ r h(x) 

smooth function of x, Xx 0, if = e ^"^ ^. 

It appears also, that class (1) contains equations being mutually equivalent with respect to 
point transformations which do not belong to this group. In particular, it is proved in [28,31,49] 
that if an equation of form (1) is invariant with respect to a Lie algebra of dimension not less than 
4 then it can be reduced by point transformations to a constant coefficient diffusion-convection 
equation (17). All such equations and corresponding transformations were found in [28,31,49]. 
Some of them were known previously [41]. For the convenience of the readers we adduce the 
results of group classification up to the extended equivalence group (22) in Appendix A. 

Up to equivalence transformations (22) the list of equations of form (1) reducible to the con- 
stant coefficient form together with corresponding transformations is exhausted by the following 
ones [28]: 

1. ut = ie^Ux)x + exux (8): i= e^^^ /{2e), x = xe^*, u = u; 

2. x~^ut = {e'^Ux)x ^ (8): i = t signx, x = \/x, u = u — \n\x\] 

3. x~^ut = {e'^Ux)x + x~'^Ux (8): i = (e^*t signa:)/2, x = e~*/x, u = u — t — ln|a;|; 

4. Ut = {\u\^Ux)x + exux (9): t = e^^^/{2e), x = xe^*, u = u; 

3m+4 , 1 

5. |x| f'+i Ut = {\u\'^Ux)x (9): t = t, X = —l/x, u = \x\ i+f 

6. jxj Ut = {\u\''ux)x + £x\x\ ''+1 n^|^^2 ^ (9): t = - e ' x = xe'\ 
u = u; 

7. x~'^ut = {u~'^Ux)x + £x~^Ux — (5): i = t, x = xe^*, u = u; 

8. e^ut = {u''^Ux)x (14): t = t, x = x, u = e^u; 

9. e^ut = iu~^Ux)x + se^Ux — (14): i = e^^/e, x = x + et, u = e^'^^^u. 
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Combining these transformations with symmetry and equivalence transformations one can 

easily obtain solutions of such "non-essentially variable coefficient" equations. For instance, 
starting from solutions of equation (8), we obtain corresponding solutions for the more compli- 
cated and interesting equation 



(76^ _|_ 1)3 

having the density / localized in space [49] : 

u = ln|ci+co(e- + 7)|, ^ = 2t(ei + 7) " T + ^"^"T^ 

Similarly one can find exact solutions of equation 

that is reducible to the fast diffusion equation [29] : 

g-(a;+£t) g.g-(a;+£t) 



U = r-rr^n-, U 



2e-^ 2e-^ 2e-^ 

U = —, r;^— rr-T, U = ttz T, U = —- 



£{x + ety ± e^^^' £cos'^{x + et)' ecosh^{x + et) 

2e-^ 2e-(^+^*) sin(2e^V£) 

U= TT. -, U = 



£:sin^(a: + et) ' cos(2e^*/^) ~ cos 2{x + et) ' 

2e-(''+-'^ siuh(2r-7c) 2(-^''-''^ siuh(2e^7s) 

u = —- 



cosh2(a; + et) - cosh(2e^V£) ' cosh2(a; + et) + cosh(2e^V£) ' 

2g-(x+£t) cosh(2e^Ve) 2e-(^+^*) sin(2e^V£) 

u = . . ^ . . ' \ , , , u = 



u 



sinh2(a; + et) — sinh(2e^*/£) ' cosh2(a; + et) — cos(2e^*/£) ' 

2g-(x+£t) sinh(2e^7£) 



cosh(2e^*/£) — cos 2{x + et) 
Using the same approach we constructed exact solutions for the following equations [29,31,49]: 

^-2x+~ie "^11^= {u~^Ux)x + U~^Ux ■ 

2cfte-T^"" _ 2teoc2e(^i-'')'^"" 



cos2ci(e-^ + co)' (1-006^1^^72' 



cie-^^^ £e-T^ 2te-T^" 
= 7—rz — 17, u= — ■ , u 



-£ + coe'=i(^""^-^*) ' e-^-et + co' (e"^ + ci)2 + cot^ ' 

Equation 

^-2x 

Ut = {u^Ux)x + u'^Ux 



4+3^ 

(e~=" + 7) i+z^ 
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has exact solutions of the form 



« = |co(e + 7) - ci|m+i, u= {cq+ +£>n |e "^ + 71 m+i , 

V + 7 ) 

Ml/ 1 \^ / 1 \ , _MaM±3) 



^ = (6t + c'i + C2e-")2(e-^ + 7)6. 

admits the following invariant exact solutions: 

/i ' |Lt cos2[ci(a; + fit + co)]' M [1 - coe^i(^+''*)]2 ' 

ci 2 p-^ 



u 



1 g-(a;+/it) 

X + //t — e'** + Co 



where / = f{x) is an arbitrary positive function has the following invariant solutions: 
u-c^,l^/jax) , ^-cos2[ci(ln/V7dx + co)]' " [1 - co(/ vTd^^F ' 

\-2 / f /7^^^-2 



= co( y v^fdx)^ 



ciij^dx)-' ^_ {jVfdxy 



u = 



C2 + Co(/ V7dz)=ie^iC2* ' ln(/ V/^^ar) - t + cq ' 

2t{f^dx)-^ 
[ln{J^dx) + ci]^ + cot^' 



admits the following solutions: 



fj, cos2 [ci (a;e'^^ + cq)] ' At(l - coe^i^e 



Cl 2jLt 

u = TT- , u — 



C2 + coe 

Aii^{xei^* + ci)2 + coe^^^* ' 

Sec [28,29,31,49] for more detail and more examples. Since the solutions of these equations 
can be reconstructed from ones presented in other sections, we tern back to the more interesting 
cases, in particular, to the equations which are "essentially variable coefficient" . 
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13 Examples of Lie reduction of variable coefficient equation 



The results of this section have been obtained in [29,31]. 

In this section we consider Lie reductions of some variable coefficient equations of form (1). 
We start from the equation 



(23) 



which admits three-dimensional Lie algebra Qi = —p~^dt, Q2 = e^'^'*{dt — d^), Q3 = dx + 
pdu- The only non-zero commutators of these operators are [Q 1,(52] = Q2- Therefore 
is a realization of the algebra 2^2. 1 [40]. All the possible inequivalent (with respect to inner 
automorphisms) one-dimensional subalgebras of A2.1 © Ai [46] are exhausted by the ones listed 
in Table 11 along with the corresponding ansatze and the reduced odes. 

Table 11. Reduced odes for (23). A 7^ 0, £ = ±1 



N 


Subalgebra 


Ansatz u — 


UJ 


Reduced ODE 


1 


Qi 




X 




2 


Q2 




t + x 


[e'^ip')' = 


3 


Qi 


(p{u)) +px 


t 


= p + p 


4 


Ql + AQ;, 


(^(^■) + Xpt 


X - Xt 


e-[Ap-(A + l)^'] = (.^^')' 


5 


Q2 -r :Ch 




X + t - -(:''' 





As a second example we consider the equation 
xPut = [u'^u^]^ ajP+^w^. 



(24) 



The invariance algebra of (24) is generated by the operators Qi = —{p + 2) ^dt, Q2 = 

is a realization of the algebra A2.1 © Ai 
too. The reduced equations for (24) are listed in table 12. 

Table 12. Reduced odes for (24). A 7^ 0, e = ±1 



N 


Subalgebra 


Ansatz u = 




Reduced ODE 


1 
2 
3 

4 

5 


Ql 
Q2 
Q3 

Ql + XQ3 

Q2 + eQi 


Lp{u}) 
Ell , . 

e^(''+2)V(c^) 
e v'(w) 


X 

xe* 
t 

xe 

xe 


{ip'^ip')' +x''^'^' = 
(^-^')' = 

,J _ (p+2)(p+2+m) m + l , p+2 , „ 
^ r m ^ 

(ST-Lp')' + (1 + Am)w''+V' 
= X{p+1)uj^ip 

(^™/)'=£(p + 2W 



At last, let us analyze in more detail equation 

eP'^'wt = (eP=^'wa=)a, + eP^'uu^, (25) 
which is invariant with respect to three-dimensional Lie symmetry algebra 

{du e-^p'dx, dx-2pdu). 

In contrast to the case of equations with four-dimensional Lie symmetry algebra wc cannot 
reduce equation (25) to a constant coefficient equation of form (1). However, it is an interesting 
feature of this equation that using a point transformation v = u + 2px it can be mapped to a 
constant coefficient reaction-convection-diffusion equation 



vt = Vxx + vvx - 2pv 



(26) 
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that does not belong to class (1). For simplification of the technical calculations we will in- 
vestigate the constant coefficient equation (26) instead of (25). The Lie symmetry algebra of 
equation (26) 

{Xi = dt, X2 = e-^P\d^ + 2pd^), Xs = d^) 

is a realization of A2.1 Ai [40]. These operators generate the following group of point trans- 
formations: 



t = t + ei, X = X + 



~2pt 



+ £3, 



V + 2e2pe 



-2pt 



A list of proper inequivalent subalgebras of the given algebra is exhausted by the following ones 

(X.+aXs), {Xs + eX2), {X2), {Xi, X3), {X2, Xs), (Xi + /3X3, X2), 

where a and /? are arbitrary constants, e = 0, ±1 [46]. 

The first three (one-dimensional) subalgebras lead to Lie reductions to ordinary differential 

equations, the fourth and sixth (two-dimensional) ones yield reductions to algebraic equations. 
Lie reductions with respect to these subalgebras are summarized in Table 13. One can easily 
check that it is impossible to construct a Lie ansatz corresponding to the subalgebra {X2, X3). 

Table 13. Lie reductions of equation (26). 



N 


Subalgebra 


Ansatz v = 




Reduced equation 


1 
2 
3 
4 
5 


(Xi + aXs) 
{X3+SX2} 

{X2) 

{Xi, X3) 
(X1+/3X3, X2) 


. , 2-pex 

(^(w) + 2'px 
C 

2px - 2p(it + C 


X — at 
t 
t 


<fi" + {'-p + a)<p' -2pip = 

= -ir <fi 

if' =G 
C = 
-2p/3 = 



Integration of equations 2-5 from Table 13 give the following invariant solutions of equa- 
tion (26): 



■u = 0, V = 2px + C, V = 



2pex + Ce^P* 



The corresponding exact invariant solutions of equation (25) have the form 

2pex + Ce^P* ^ 
u = -2px, u = C, u= — ^2pt _)_ g ^P^^ 

where C is an arbitrary constant. 

Ansatzes 4.2 and 4.3 give a hint for a possible form 



(27) 



V = ip{t)x + V'(i) 

of nonlinear separation of variables for construction of exact solutions of equation (26) . Substi- 
tution of the ansatz to equation (26) leads to antireduction: 

(p' = (p'^ — 2pip, if)' = (ftp — 2p^. 

Solving the above system of ODEs for (p and ip we obtain exactly the solutions of equations 4.2 
and 4.3. 
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Note 1. Using the point transformation i = e ^P^, x = x, v = e^P^v equation (26) can be 
mapped to a variable coefficient Burgers equation = —2pt~^Vxx — ^pvv^ studied in [37]. 

Note 2. The well-known Cole-Hopf transformation v = 2wx/w reduces equation (26) to the 
famous constant coefficient reaction-diffusion equation with weak nonlinearity 

u>t = Wxx — 2pw \n\w\. 

After application of the Cole-Hopf transformation to the list of known exact solutions (see, 
e.g., [48]) of the equation with weak nonlinearity we obtain exactly solutions (27) of equa- 
tion (26). 



14 Exact solutions of sZ(2, ]R)-invariant equation 

Analyzing the results of group classification of diffusion-convection equations, we can observe a 
number of G~-inequivalent equations (1) which are invariant with respect to different realizations 
of the algebra s/(2,]R). The set of such equations is practically exhausted by the well-known 
("constant coefficient") Burgers and «~^/^-diffusion equations and by the equations which are 
equivalent to them with respect to additional transformations. This set is supplemented by the 
unique essentially variable coefficient equation [28, 31] 

x^ut = {u-^/^Ua,)a, + x^u^, (28) 

s/(2, ]R)-invariance of (28) is directly connected with the fact that h is not constant. The 
corresponding realization 

Pt = dt, D = 2tdt + 2xdx - 5udu, U = t'^dt + {2tx + x'^)dx - 5{t + x)udu. 

of the algebra sZ(2,R) is quite different from ones of cases of Burgers and -u^^/^-diffusion equa- 
tions and is the maximal Lie invariance algebra of equation (28). It was a reason to study 
equation (28) from the symmetry point of view in detail in [29]. These operators generate the 
following one-parameter groups of point transformations: 

Pt'. t = t + e, X = x, u = u; 
D: t = e^t, X = e^x, v = e^^v; 

TT -~ iy ^ i/ ~\~ X t / / \ \ (S 

n: t= , x=- 7 7 — z , u = (1 — £(t + x))u. 

The complete Lie invariance group is generated by both the above continuous transforma- 
tions and the discrete transformation of changing of sign in the triple {t,x,u). The transforma- 
tions from G™^^ can be used for construction of new solutions from known ones. 

A list of proper G^'^^-inequivalent subalgebras of is exhausted by the algebras {Pt), 

(D), {Pt + n), {Pt, D). Reduction of (28) with respect to these subalgebras and application of 
the invariance transformations lead to the following set of G°^^-inequivalent Lie invariant exact 
solutions (below 5 € {0, 1}): 



24 



In [29] it was proposed to use functional separation of variables 
/ 6 




u=iJ2^'{t)x'] . (29) 

\i=0 

to obtain solutions of equation (28). The set of all solutions of the form (29) is closed with 
respect to transformations from G™^^ and is exhausted, up to translations with respect to t and 
scale transformations, by the above solutions u = 5 and u = 5(t + x)~^ and the solutions given 
by the generalized ansatz 

U = (2x3 + ip\t)x'^ + i/{t)x^ + /(i)^6^-5/6 

and the corresponding reduced system 

^t = ^^'-\{^')\ ^l = lS^'-\^'^\ ^? = -^(^^)2 + 2//. (31) 

System (31) can be reduced to the single third-order ordinary differential equation on the func- 
tion ip'^: 

634, + ^^n^tf + 126/4 + 192(4) Vt + 16(4)^ = 

having particular solutions 4 = 0, 4 = C/t, where C G {0,3/4,21/4,6}, that lead to Lie 
invariant solutions of (28). 
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A Group classification of diffusion— convection equations 



Table 14. Case of VA{u) (gauge g = 1) 



N 


B{u) 


/(^■) 


;i(a;) 


Basis of A'"^" 


1 


V 


V 


V 


dt 


2a 


V 




1 


dt, ptdt + dx 


2a' 


V 


\x\P 




dt, {p + 2)tdt + xdx 


2b 


1 




e" + /3 


dt, e-\dt-d^) 


2c 


1 




xla;|P + /3a;-i 


dt, e-(p+2)*(c>t-xa,) 


3 


1 




x""' In \x\ 


dt, e~^xdx 


4 





1 


1 


dt, dx, 2tdt + xdx 



Here p € {0, 1} mod in case 2; p / — 2 in case 2c; (3 £ {0, ±1} in case 2b. 
Additional equivalence transformations: 

1. 2a(p = 0, B = 1) ^ 2a(p = 0, B = 0): i = t, x = x + t, u = u; 

1'. 2a' (p = -2, B = I) ^ 2a'(p= -2, B = 0): i = t, x = xe\ u = u; 

2. 2b 2a{B = P,p=l): t = e* , x = x + 1, u = u; 

3. 2c(p^ -2) -» 2'a{p^-2): i= {e'-P+^'>^ - l)/(p + 2), x = xe\ u = u. 



Table 15. Case of A{u) = e" 



N 


B{u) 




9{x) 


h{x) 


Basis of A""^" 


1 





V 


1 


1 


dt, 


tdt - du 


2 


e"" 




1 


\xy 


dt, 


{pn - pu - 2v - qn + n)tdt + (m - '^)xdx + {q+ l)du 


2* 






1 


ee" 


dt, 


{pn -pv - i^)tdt + i')dx + du 


3 




^px'^+qx 






dt. 


{2p + q)tdt + dx- 2pdu 


4 




1 


1 


1 


dt. 


dx, (p — 2v)tdt + ((/i — y)x + vxt)dx + du 


5 




1 


1 


1 


dt. 


dx, tdt + {x- t)dx + du 


6a 





/Ha;) 


1 


1 


dt, 


tdt — du, atdt + {Px^ + 7ix + 7o)9a! + Pxdu 


6b 


1 


\x\^ 


1 


ex\x\^ 


dt. 


xdx + b + 2)du, e-''^^+^^\dt - sxdx) 


6b* 


1 




1 


ee" 


dt. 


dx+du,e-'\dt-edx) 


6c 


1 




1 




dt. 


xdx, tdt — stxdx — du 


7a 





1 


1 


1 


dt, 


tdt — du, 2tdt + xdx, dx 


7b 


1 


1 


1 


1 


dt. 


dx, tdt - td,r - d,„ 2tdt + (x - t)dx 


7c 


1 


1 


1 


ex 


dt. 


xdx + 2du, e-^'dx, e-''''{dt - exdx) 


7d 





x-^ 


1 


1 


dt, 


tdt — du, xdx — du, x^dx + xdu 


7e 


1 




1 


a;-2 


dt, 


xdx — du, e*{dt — xdx), e*{x'^dx + xdu) 



Here (p, I/) G {(0, 1), (1, v)}, v ^ n in cases 2, 2* and 4; /x = 1 and v ^ 1 in the other cases; q 7^ —1 in 
case 2* (otherwise it is subcase of the case 1.2a'); e = ±1 in cases 2, 6b-6c and 7e; p ^ {—3, —2, 0} in case 6b; 
a, /3, 71 , 70 = const and 

-3/3x — 271 + a 



f (x) = exp 



■ dx 



l3x^ + jix + 70 

Case 2{q = —1) is a subcase of case 1.2a'. Additional equivalence transformations: 

1. 4:{x 7^ 0) ^ 4(x = 0): i = t, X = X + xt, u = u; 

2. 6b -> 6a (/3 = 70 = 0, q;= (p + 2)7i): t= (e^(P+'^^* - l)/(e(p + 2)), x = xe^\ u = u; 

3. 6b* — » 6a (/3 = 71 = 0, a = 70): t = e^*/e, x = x + st, u = u; 

4. 6c ^ 6a (/3 = 70 = a = 0): i = t, x = xe^*, u = u; 

5. 7b— >7a: i — t, x — x + t, u = u; 

6. 7c^7a: i = e^=V(2e), x = xe^\ u = u; 

7. 7d— »7a: t = tsigna;, x = l/x, fi = u — ln|a;|; 

8. 7e— »-7a: i = {e^*tsignx)/2, x = e~*/x, u = u - t - ln\x\. 
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Table 16. Case of A{u) = \u\' 
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c 
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8 


V 




1 


1 


1 


dt, dx, 


Q 


w 

V 


In U/l 
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1 

J. 


1 

L 


1 

J. 


Ot, <Jx, t^iut yH'X — iJOx uuu 
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It 


6 
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C 


e 


Ot, e Ox, Ox — ^pOu 


il 


U 


In L.I 
111 \U\ 








Ut, e Ox, Ox — 2pUOu 


12a 


V 





fix) 


1 


1 


dt, fJ.tdt — udu, 

aiOt + [(n + ^)px + 7ia: + 'yojOx + pxuOu 


1 9K 
iZD 


w 
V 


1 
1 


|x| 


1 
1 




Ot, uxOx + \JP + 'ijuOu, e (Ot — sxOxj 


Izu 




1 

i 


a: 

e 


1 
i 


X 

ee 


Ot , l-lOx + uOu , e {Ot — £Ox ) 


IZC 




1 

± 




1 

± 






13 


—6/5 


1 




1 




t^Q^ + ^2tx + x^)dx — 5{t + x)udu 


14a 


/-4/3 





1 


1 


1 


dt, fJ,tdt — udu, dx, 2tdt + xdx 


14b 


/-4/3 


1 


1 


1 


1 


dt, )J,tdt — ^J,tdx — udu, dx, 2tdt + (x — t)dx 


14c 


7^-4/3 


1 


1 


1 


ex 


dt, nxdx + 2udu, e-^'dx, e-''^\dt - exdx) 


14d 


^ -4/3, -1 





|a;| ''+1 


1 


1 


dt, i^tdt - udu, in + 2)tdt - (/i + Vjxdx, 








|a;| f'+i 




£x\x\ '■+1 


(/i + l)x^dx + xudu 


14e 


-4/3, -1 


1 


1 


dt, m(m + '^)xdx — + 2)udu, 

e ''+1 \at—exOx),e ((/Lt+ijx Ox-yxuOu) 


14f 


-1 







1 


1 


dt, tdt + udu, dx — udu, 2tdt + xdx — xudu 


14g 


-1 

— i 


-1 

1 


e 


-1 
i 


£6 


Ot, Ox — UOu, (X + et — 2)Ox — [X + etjuOu, 

e-''{dt-edx) 


1 /I Vi 


— Z 


1 
i 


X 


1 
i 




Ot, xOx, ztot — zstxOx + UOu, 
e^*{x^dx — xudu) 


ioa 


-4/3 


U 


1 


1 
i 


1 


Ot, 'aot + ouOu, Ox, Mot + xOx, 
x^dx — .ixudu 


iOD 




1 
1 


1 
1 


1 
1 


1 
1 


Ot, ^lOt -r ^XOx — oUOu, ^lOt -]- \X — tjOx, 

dx, {x + t^dx — 3(x + t)uBu 


15c 


-4/3 


1 


1 


1 


£X 


dt, 2xdx — Sudu, e~^*dx, 

e-^^\dt - exdx), e^'ix'^dx - Zxudu) 
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u 


1 


1 


1 


dt, dx, tdx - du, 2tdt + xdx - udu, 
t^dt + txdx — {tu + x)du 



Here v ^ e = ±1; q ^ —1 in case 2* (otherwise it is subcase of the case 1.2a'); p ^ —2, — (3/i + 4l)/{ix + 1) in 
case 12c; a, /3, 71, 70 = const, and 



f{x) = exp 



-(3m + 4)/?2:-3 



dx 



{p + l)/3x^ + X 
Additional equivalence transformations; 
1. 8(X7^0) -» 8(x = 0): i=t,x = x + >et,u = u; 



f (x) = exp 



-(3At + 4)/3a;-27i+Q 
{fl + l)Px^ + 71a; + 7o 



dx 
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2. 12b -> 12a (/? = 70 = a = 0), 14e -» 14a: i= {e^P+^^* - l)/(p + 2),x = xe*, u = u; 

3. 12b* ^ 12a (/3 = 71 = 0, a = 7o): t = (e^(P+2)* - l)/(£(p + 2)), x = xe^\ u = u; 

4. 12c 12a (/3 = 70 = a = 0), 14h ^ 14a: t = t, x = xe^\ u = u; 

5. 14b — > 14a, 15b — > 15a: i = t, x = x — t, u — u; 

6. 14c -» 14a, 15c -» 15a: f = e^^7(2£), x = xe^*, ■« = «; 

7. 14d -» 14a: f = t, x = -1/a;, u = \x\~t^u; 

8. 14f — > 14a(// = -1): f = t, a = x, w = e^'w. 

9. 14g -» 14a(/i = -1): i = e^*/e, x = x + et, u = e^+^*w; 
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